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Among the phases of SiC>2 are a- and /3-cristobalites, which have a long and somewhat controver- 
sial history of proposed structural assignments and phase-transition mechanisms. Recently, Zhang 
and Scott found new indications that the higher-temperature p phase has space group lA2d and, by 
assuming a group-subgroup relationship between phases, they argued that the lower-temperature 
a phase should have lower symmetry than that of the widely-accepted P4i2i2 space group. With 
this motivation, we use first-principles calculations to investigate the energy, structure, and local 
stability of P4i2i2 and J42d structures. We also compute the frequencies of the zone-center phonon 
modes in both structures, as well as certain zone-boundary modes in the I42d structure, and com- 
pare with experiment. We then argue that the various P4i2i2 and I42d enantiomorphs can be 
grouped into three clusters, each of which is identified with a three-dimensional manifold of struc- 
tures of P2i2i2i symmetry in which the P4i2i2 and I42d appear as higher-symmetry special cases. 
We find that there are relatively high energy barriers between manifolds, but low barriers within a 
manifold. Exploring the energy landscape within one of these manifolds, we find a minimal-energy 
path connecting P4i2i2 and I42d structures with a surprisingly low barrier of ~5meV per formula 
unit. Possible implications for the phase-transition mechanism are discussed. 

PACS numbers: 61.66.Fn, 63.20.dk, 64.60.Ej 



I. INTRODUCTION 

The fact that SiC>2 can exist in numerous crystalline 
and amorphous forms, and its status as one of the most 
prevalent compounds on earth, has stimulated a long his- 
tory of experimental and theoretical investigation. Here 
we focus on the a ("low") and (3 ("high") cristobalite 
phases, which are stable near the melting temperature 
and are metastable at room temperature. 

The structure of the higher-temperature f3 phase has 
a history of controversy. Early indications of a cubic 
structure with 180° bond angles (space group FdSm) 1 
were challenged by others 2,3 who hypothesized that the 
true /3-phase structure has lower symmetry and that 
the apparent cubic structure arises from averaging over 
spatial domains or dynamical fluctuations. In partic- 
ular, Wright and Leadbetter 3 argued for a tetragonal 
structure belonging to space group I42d (D^). While 
some subsequent works have provided support for this 
identification, 4,5,6 other authors maintain that it is better 
to describe the (3 phase as a dynamically disordered one 
having overall Fd3m symmetry but with a large popula- 
tion of rigid-unit-mode (RUM) fluctuations. 7,8 To some 
degree, the argument may be semantic; if the fluctua- 
tions have strong short-range correlations in space and 
time, it is difficult to distinguish this picture from one of 
dynamic domains of a lower-symmetry structure. 9 Thus, 
for example, either picture may be able to explain the 
fact that there are two more first-order lines in the Ra- 
man and infrared spectra than would be expected from 
Fd3m symmetry, 6 and the question of which description 
is "correct" might depend on the time and length scales 
of the experimental probes in question. 

In contrast, the assignment of the a-cristobalite phase 
to the tetragonal P4i2i2 (Df) space group 10 has un- 



til recently been noncontrovcrsial. However, based on 
a reexamination of Raman and infrared vibrational spec- 
troscopies, Zhang and Scott 6 have recently raised new 
questions about the identity of the a phase. By using 
Raman spectroscopy to study small single crystals of (3- 
cristobalitc, these authors argued that the f3 structure 
must be D 2 d, not cubic, and assuming a group-subgroup 
relationship for the /3-to-a transition, concluded that the 
a phase should have some lower symmetry such as D2 in- 
stead of D4. The apparent D4 symmetry of a-cristobalite 
could result from spatial or dynamic averaging over D2 
domains, in analogy to what had been proposed for the 
[3 phase. To support their assumption that a group- 
subgroup relationship should hold, Zhang and Scott also 
pointed to the temperature dependence of the optical 
phonon frequencies near the phase transition as being 
inconsistent with a reconstructive phase transition 11 and 
as suggesting a nearly second-order behavior, although 
arguing in the opposite direction are the facts that the 
latent heat and volume change at the transition are quite 
substantial. 12 

In their paper, Zhang and Scott 6 reexamined ear- 
lier Raman and infrared spectroscopic measurements not 
only on the a- and /3-cristobalite SiC>2, 13 but also on a 
and f3 AIPO4 (Ref. 14) and a BP0 4 (Ref. 15) cristo- 
balites. Note that the replacement of Si atoms by Al and 
P (or B and P) atoms immediately reduces the symme- 
try according to P4A2 (Dj) -> C222x (Df) for the a 
phase 16 and H2d (Z? 12 ) -> 14 (Sf ) for the phase. Also 
of possible relevance is the pressure-induced phase tran- 
sition from cv-Si02 to a high-pressure monoclinic cristo- 
balite phase. 17 The relationship of these other cristo- 
balites to the a and (3 phases of SiC>2 is an interesting 
avenue for future exploration, but falls outside the scope 
of the present work. 

First-principles calculations of the structural and lat- 
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tice dynamical properties of SiC>2 have a long and pro- 
ductive history. 5 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 ' 30 ' 31 While 
quite a few of these works specifically address the a- 
cristobalite structure, 18 ' 19 ' 20 ' 24 ' 25 ' 27 ' 31 questions about its 
stability and about possible pathways from the a to the 
(3 phase have not been fully explored. 

In the present work, we have carried out first-principles 
calculations for candidate a and (3 cristobalitc structures 
in the framework of density-functional theory (DFT) in 
order to check the stability of both phases and to ex- 
plore the energy landscape connecting them. We have 
also calculated phonon frequencies and infrared activ- 
ities for both a and f3 phases, and explored how the 
phonon modes in the different phases are related to each 
other and to those of the high-symmetry cubic phase. 
Our calculations are effectively zero-temperature ones, 
and thus cannot properly treat the temperature-induced 
a-(3 cristobalite phase transition. Nevertheless we hope 
that the information obtained from these calculations can 
eventually be built into a realistic statistical-mechanical 
theory that correctly takes the RUM fluctuations into 
account in its description of the a and (3 phases at ex- 
perimentally relevant temperatures. 

The paper is organized as follows. In Sec. II we give 
a brief review of a and (3 cristobalitc structures and de- 
scribe the methods used in the calculations. Then, in 
Sec. Ill, we present the results of our calculations of struc- 
tural and lattice vibrational properties of the two phases 
and of the energy landscape connecting them. We dis- 
cuss those results in Sec. IV. Finally, we summarize the 
work in Section V. 



II. PRELIMINARIES 

A. Cristobalite structures 

In order to describe the structures of the SiC>2 a and 
(3 cristobalitc phases, it is easiest to start by considering 
the "ideal cristobalite" structure, which is constructed 
by placing Si atoms in a diamond structure with oxy- 
gen atoms located midway between each pair of nearest- 
neighbor Si atoms. This structure has the space group 
Fd3m (0^) and has two formula units per primitive unit 
cell. Each Si atom with its four surrounding O atoms 
forms a tetrahedron, so the whole structure can be vi- 
sualized as a network of tetrahedra connected at their 
apices. 

The generally accepted structure of a-cristobalite is 
arrived at by starting from the ideal structure and mak- 
ing nearly rigid rotations of the tetrahedra about [100] 
and [010] axes (in the original diamond cubic frame), 
leading to a tetragonal structure with its axis along z. 
This is illustrated in Fig. 1(a), but in the conventional 
tetragonal frame, related to the original cubic frame by 
a 45° rotation about z. The tetrahedral rotations are 
also accompanied by small strains and tetrahedral trans- 
lations needed to keep the apices coincident, as would be 
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FIG. 1: (Color online) Projection on x-y plane of the (a) 
a and (b) (3 structures, proposed as candidates for a and /3 
cristobalite phases respectively. Darker shading is used to 
represent more distant tetrahedra so that the spiral structure 
of the connected tetrahedrons becomes evident; double ver- 
tical lines indicate that the adjoining tetrahedra are actually 
disconnected because they are separated in the z-direction. 



expected from enforcement of the rigid-unit constraints. 
The space group of the structure is P4i2i2 (Df ), and 
since the four rotations shown in Fig. 1(a) are all differ- 
ent, the number of formula units per primitive unit cell 
is now increased to four. 

As mentioned earlier, diffraction experiments on the 
/3-cristobalite phase tend to give inconclusive results be- 
cause of spatial and dynamical averaging. Nevertheless, 
based on the comparison of structure factors predicted by 
various disorder models and the ones obtained in their x- 
ray diffraction experiments, Wright and Leadbetter 3 con- 
cluded that the /3-cristobalite has local I42d (D^) space- 
group symmetry. Their proposed structure can also be 
constructed from the ideal structure, but this time by ro- 
tating all the tetrahedra around the z axis, yielding the 
structure shown in 1(b). The number of formula units per 
primitive unit cell remains at two as in the ideal struc- 
ture (although it can alternatively be described, as in 
Fig. 1(b), by a doubled conventional cell containing four 
formula units). Again, the structure is highly consistent 
with the rigid-unit constraints. 

Because we do not want to presuppose an identifica- 
tion of a particular experimentally observed phase with a 
particular crystal structure, we henceforth adopt a nota- 
tion in which the phases are identified by labels "a" and 
"/?" without tildes, whereas the putative crystal struc- 
tures shown in 1(a) and (b) will be referred to as "a" 
and "/?" structures, respectively. Our working hypothe- 
sis is that the a and (3 phases have microscopic crystal 
structures of type a and (3 respectively, but we adhere 
to a distinction in the notation in order to discriminate 
clearly between the specified structures used in our calcu- 
lations and the hypothetical identification of these with 
experimental phases. 
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B. Computational methods 

The calculations were carried out using the ABINIT 
implementation 32 of density-functional theory with 
Perdew-Burke-Ernzcrhof 33 version of the generalized gra- 
dient approximation (GGA) for electron exchange and 
correlation. Since it is the smallest unit cell that con- 
tains both a and structures, all calculations were 
performed on the four-formula-unit computational cell 
shown in Fig. 1, even though the primitive cell is smaller 
in the structure. The Brillouin zone was sampled by 
a 4 x 4 x 4 Monkhorst-Pack grid. 34 Structural proper- 
ties were computed using projector augmented-wave 35 
potentials converted from ultrasoft pseudopotentials 36 
with a plane-wave cutoff of 22 Ha unless otherwise 
specified, while phonon frequencies, eigenvectors, and 
Born charges were computed 37 using norm-conserving 
Trouillcr-Martins pseudopotentials 38 at an energy cutoff 
of 50 Ha (after repeating the structural relaxation using 
these potentials) . The acoustic sum rule was imposed on 
the force constants, and charge neutrality was imposed 
on the Born charges. Throughout the paper, the sym- 
metry analysis associated with crystal space groups has 
been carried out using the Bilbao package. 39,40 



III. RESULTS 

A. Structural properties of a and structures 

We started our calculations by considering the ideal 
cubic structure and relaxing its volume, obtaining 
a c =7.444A for the lattice constant of its eight-formula- 
unit cubic cell. Then, working in the frame of the four- 
formula-unit tetragonal cell, we found the phonon fre- 
quencies at the r point of its Brillouin zone, correspond- 
ing to phonons at the T point and one X point [namely 
(27r/a c )(001) or cquivalcntly (27r/a c )(110) in the cubic 
frame] of the primitive two-formula-unit fee cell. For the 
"ideal structure" of space group Fd3m, the symmetry 
decomposition of these phonons into irreducible repre- 
sentations is 



T (ideal) = 1A 2U © 1E U 
X (ideal) = 3Xi ffi 1X 2 



B 2T lu ffi 1T 2U 
2X 3 © 3X 4 . 



1T 2 



(1) 
(2) 



(The translational Ti u mode has been omitted.) The E u 
mode and all X modes are doubly degenerate; the Ti u , 
T 2u and T 2g modes are triply degenerate; and A 2u is 
non-degenerate . 

We found that the triply-degenerate T 2u mode at T is 
unstable with an imaginary frequency of «83cm _1 . All 
other optical phonons have real frequencies, the lowest 
being at 250 cm -1 . Furthermore, one of the doubly- 
degenerate (X4) modes is unstable with a frequency of 
ibScm" 1 . We thus conclude that the ideal cristobalitc 
structure is unstable with respect to these distortions. 
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FIG. 2: Ground state energy per formula unit (eV) vs. vol- 
ume per formula unit (A 3 ) for a, (3, and cubic cristobalite 
structures of Si02- 



Next we imposed distortions corresponding to these 
unstable modes and did a full relaxation of the struc- 
ture subject to the symmetry constraints of the resulting 
space group. The unstable (z53 cm -1 ) mode at X leads to 
the space group P4i2i2 (or P4a2i2) which corresponds 
to ev-cristobalite, while the (i83cm _1 ) mode at T takes 
us to the space group I42d of /3-cristobalite. The en- 
ergy of the relaxed ground state as a function of volume 
per formula unit is shown for both cases in Fig. 2, with 
the energy of the cubic phase also shown for reference. 
The corresponding structural parameters at the energy 
minimum are given in Table I. 

From Fig. 2 it is clear that the a and structures in- 
deed have lower energies than the ideal cristobalite when 
the volume becomes smaller than some critical volume 
V - 55 A 3 . (Ab ove this volume, the imposed distortions 
disappear during relaxation and the structure returns to 



TABLE I: Lattice constants (in A) and Wyckoff structural 
parameters for a (P4i2i2) and (3 (742d) cristobalite Si02. 





Present 
theory 


Previous 
theory" 


Expt. 6 


a-cristobalite 








a 


5.0730 


5.1190 


4.9570 


c 


7.0852 


7.1683 


6.8903 


Si{u) 


0.3001 


0.2869 


0.3047 


O(x) 


0.2384 


0.2439 


0.2381 


0(2/) 


0.1081 


0.0777 


0.1109 




0.1819 


0.1657 


0.1826 


/3-cristobalite 








a 


7.1050 


7.226 


7.13L 


c 


7.4061 


7.331 


7.13L 


O(x) 


0.1051 


0.0896 


0.079 



°Ref. [24]. 

6 Refs. [3] and [10]. 

c Experiment sees average cubic structure. 
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the ideal one.) We find that both the a and (3 structures 
have a quite similar dependence of energy on volume. Ac- 
cording to our calculation, the relaxed (3 structure has a 
slightly lower energy than that of the a structure (12 meV 
per formula unit). This appears to be in conflict with 
the experimental situation, since the a phase is exper- 
imentally more stable at lower temperatures. However, 
when we repeated our calculations using a local-density 
approximation (LDA) exchange-correlation functional, 41 
the (3 structure was found to be lower by 1 meV per 
formula unit. We thus conclude that the small energy 
difference between the two nearly-degenerate structures 
is a quantity that is too delicate to be reliably obtained 
by our DFT calculations. A similar discrepancy between 
the results from LDA and GGA functionals was found in 
Rcf. [24]. 

We also analyzed what happens to the bond lengths 
and angles in the a and (3 structures as a function of vol- 
ume. The details of this analysis are deferred to App. A. 
Briefly, for V < Vq, the O-Si-0 bond angles and Si-0 
bond lengths inside the tetrahedra are found to remain 
almost constant, while the Si-O-Si bond angles change 
by ^35°. The details are shown to be very close to the 
predictions of a picture of tilting of perfectly rigid tetra- 
hedra. The fact that the three phases are indistinguish- 
able for V > Vo is also easily explained, since the tilts 
of rigid tetrahedra can only decrease the volume of the 
ideal structure. Thus, for V > Vq the tetrahedra cannot 
remain rigid and the Si-0 bond length must increase, 
and only when V becomes smaller than some volume Vo 
will one of the RUMs condense in the structure in order 
to maintain the preferred Si-0 bond lengths. 



B. Phonons 

1. Phonons at T in a cristobalite 

We next repeated the calculation of the phonon fre- 
quencies for the fully relaxed a and (3 cristobalite struc- 
tures. We did this in order to compare with experimental 
measurements and also to check the stability of the struc- 
tures and to investigate, at least at harmonic order, the 
nature of the energy landscape around these structures. 
An analysis extending beyond the harmonic approxima- 
tion will be presented in Sec. IIIC 2. 

The decomposition of the optical Y phonons into irre- 
ducible representations for the a structure in space group 
P4i2i2 is 

T (5) = 4Aj © 4A 2 © 5Bi © 4B 2 © 8E. (3) 

(The translational A 2 and E zero modes have been omit- 
ted.) Only the E modes are doubly degenerate; all others 
are non-degenerate. 

Tables II and III present the phonon frequencies at the 
r point for the fully relaxed ci-cristobalitc structure. All 
phonon frequencies are positive, although some appear to 



TABLE II: Infrared-active phonon modes at L in a- 
cristobalite (P4i2i2). (E modes are also Raman-active.) For 
A2 modes, ulo refers to a phonon with q = z, while for E 
modes ojlo refers to q lying in x-y plane. 



Irrep 


(l)rT^f^ ( PTTI ^ ^ 




!i)t 1-4 ( ( 

^LO \* 


E 


127 


0.05 


128 


E 


259 


0.04 


260 


A 

A 2 


ZOO 


n on 


zyo 


E 


357 


0.18 


360 


E 


440 


0.74 


507 


A 2 


462 


0.67 


515 


E 


584 


0.23 


591 


A 2 


751 


0.52 


764 


E 


752 


0.02 


752 


A 2 


1050 


1.52 


1201 


E 


1170 


0.17 


1165 


E 


1048 


1.55 


1208 



be rather low in frequency. For the infrared (IR) active 
modes shown in Table II, the transverse mode frequen- 
cies were computed initially, and their mode dynamical 
charges were also computed using 

Za, Q = ^2 ~lj^^P Z ta0 ( 4 ) 

iff v 1 

where £i \p is an eigenvector of the dynamical matrix, 
Z* a p is the Born atomic charge tensor, and Mi is the 
mass of the i-th atom in amu. The norms of the mode- 
charge vectors Z A = E„(Zj a ) 2 ] 1 '' 2 are also given in 
the Table. The longitudinal dynamical matrix was then 
constructed and diagonalized using standard methods, 37 
and the resulting LO mode frequencies are presented in 
the last column of Table II. It can be seen that there 
are wide variations in the mode dynamical charges, and 
consequently, large variations in the LO-TO splittings. 

TABLE III: Raman-only phonon modes at V in a-cristobalite 



(P4i2i 


2). 










Irrep 


iv (cm -1 ) 


Irrep 


u! (cm 1 ) 


Irrep 


u) (cm 1 ) 


Bi 


29 


Bi 


358 


Ai 


1046 


Bi 


103 


Ai 


378 


Bi 


1049 


Ai 


197 


B 2 


410 


B 2 


1109 


B, 


275 


Bi 


745 






Ai 


350 


B 2 


750 







2. Phonons at T and M in f3 -cristobalite 

Similar calculations of phonon frequencies were also 
carried out for the (3 structure proposed by Wright and 
Lcadbetter 3 for /3-cristobalitc. Since the primitive cell 
of the (3 structure contains only two formula units while 
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TABLE IV: T phonons in /3-cristobalite (742d) . In italics we 
show for each phonon in /3-cristobalite a closest phonon in 
a-cristobalite (P4i2i2). 



Phonon in /3 structure 


Closest in a structure 


Irrep uj (cm 1 ) 


Za 


Irrep 


uj (cm 1 ) 


Infrared and Raman 








E 126 


0.06 


E 


127 


B 2 425 


0.79 


A. 


462 


E 444 


0.72 


E 


440 


E 748 


0.51 


E 


752 


B 2 1038 


1.56 


A j 


1050 


E 1047 


1.52 


E 


1048 


Raman only 








Ai 289 




Bt 


29 


Bi 406 




A, 


350 


Bi 737 




lh 


745 


Inactive 








A 2 357 




Bg 


410 


A 2 1097 




B_> 


1109 



the a structure contains four, it should be kept in mind 
that the V point of the a structure maps not only into 
the r point of the (3 structure, but also into a second 
point that would be denoted asX= (27r/a c )(110) in the 
original fee frame, or (27r/a)(100) (where a ~ a c /^/2) in 
the rotated frame of Fig. 1; we shall refer to this as the M 
point in accordance with the conventional labeling of the 
body-centered-tetragonal (bet) primitive cell in the latter 
frame. The decompositions of the T and M phonons into 
irreducible representations for the f3 structure in space 
group IA2d are 

T (J\ = lAi © 2A 2 © 2Bi © 2B 2 © 4E, (5) 

M (p\ = 2MiM 2 © 3M 3 M 4 © 4M 5 . (6) 

(Translational B 2 and E zero modes have been omitted.) 
All M-point modes and T-point E modes are doubly de- 



TABLE V: M phonons in /3-cristobalite (742d). In italics we 
show for each phonon in /3-cristobalite a closest phonons in 
a-cristobalite (P4i2i2). 

Phonon in /3 structure Closest in a structure 

Irrep uj (cm -1 ) Irrep uj (cm -1 ) Irrep uj (cm -1 ) 



M3M4 35 Bi 103 Ai 197 

M 5 281 E 259 E 357 

MiM 2 316 B 2 275 A 2 285 

M3M4 336 B t 358 A: 378 

M 5 372 E 259 E 357 

M 5 586 E 584 

MiM 2 780 B a 750 A 2 751 

M3M4 1045 A: 1046 Bi 1049 

M 5 1162 E 1170 



TABLE VI: Relations between unstable phonons in "ideal 
structure" and phonons in a and (3 structures. 





Ideal 


a-cristobalite 


/3-cristobalite 




cm -1 


cm 1 


Irrep 


cm 1 


Irrep 


r 


i83 


29 


Bx 


289 


Ai 


r 


i83 


127 


E 


126 


E 


r 


i83 


127 


E 


126 


E 


X 


i53 


197 


Ai 


35 


M3M4 


X 


i53 


103 


Bi 


35 


M3M4 



generate, while other modes are non-degenerate. 

In Tables IV and V we present the our results for 
the L-point and M-point phonon modes, respectively, of 
/3-cristobalite. The frequencies given for the IR-active 
modes at T are the transverse ones only. The tables 
also show the correspondences between the T modes in 
the a-cristobalite structure and the T and M modes in 
the /3-cristobalite structure, as determined by comparing 
phonon eigenvectors. 



3. Relation to unstable phonons in the cubic phase 

The triply degenerate T-point mode of the cubic struc- 
ture having imaginary frequency i83cm _1 , which con- 
densed to form the (3 structure, now has positive frequen- 
cies of 289 cm -1 for the non-degenerate Aj mode and 
126 cm -1 for the E doublet in the f3 structure. This same 
triplet corresponds, in the a structure, to the lowest- 
frequency phonon of frequency 29 cm -1 , which has sym- 
metry Bi, and to an E doublet at 127 cm -1 having al- 
most the same frequency as in the (3 structure. The 
doubly-degenerate unstable mode of the cubic structure 
at i53cm -1 , which condensed to form the a structure, 
now appears in the a structure at frequencies 197 cm -1 
and 103 cm -1 with symmetries Ai and Bi, respectively. 
In the (3 structure, on the other hand, the same doublet 
appears as the lowest-frequency phonon in that structure, 
namely the doublet at 35 cm -1 with symmetry M3M4. 
These relations between the unstable phonons in the 
"ideal structure" and the phonons in a and f3 structures 
are shown in Table VI. 



4- Comparison with experiment for a-cristobalite 

In view of the recent questions that have been posed 
about the identity of the a-cristobalite phase, 6 we have 
carried out a more detailed analysis of the phonons in 
the a structure. In particular, we have calculated the LO 
frequencies of the T-point phonons in a-cristobalite as a 
function of the angle at which the limit g — s- is taken. 
It turns out that the labels A 2 and E are not well-defined 
at arbitrary q because of mixing between modes of these 
symmetries. Moreover, it can happen that if one starts 
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TABLE VII: Left: Computed mode frequencies and irreps for a-cristobalite, with direction of dynamical polarization in 
parentheses for IR-active modes. The modes that are adiabatically connected as q is rotated from x to z appear on the same 
line. All modes other than A2 modes are Raman-active. Right: Tentative assignments to measured mode frequencies in powder 
samples. 
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a Not pure LO at all q. b LO-TO splitting is negligible. c Part of structured peak. 
d Inactive in /3 phase. In (3 structure the A 2 component also becomes Raman active. 
f Corresponds to RUM mode in ideal cristobalite. s Observation from [13]. h Observation 
of SOcm^ 1 mode is from [42], all others from [43]. 1 Empirical assignments from [13]. 



with an E mode at q \\ z and follows the branch as q is 
rotated, one arrives at an A2 mode when q lies in the x- 
y plane, or vice versa. Experiments have typically been 
done on powder samples, so that one should in principle 
average the phonon spectrum over all possible directions 
for q — > 0. Moreover, some phonon modes with E sym- 
metry have a very small LO-TO splitting, so they would 
most likely appear in experiment as a single line. 

For all these reasons, a direct comparison of experi- 
mental data with our results as presented in Tables II and 
III is problematic. Nevertheless, we attempt such a com- 



parison in Table VII. Despite the difficulties, the agree- 
ment with experimental data is rather good, with a few 
exceptions that will be discussed shortly. We generally 
underestimate the experimental frequencies by ~20 cm -1 
for lower frequency phonons and by ~35cm _1 for higher 
frequencies, but otherwise our results reproduce the ex- 
perimental pattern of frequencies, and the irrep assign- 
ments arc also consistent with those obtained from em- 
pirical models. 13 Furthermore, the identification of the 
modes that are not expected to be active in the (3 phase 
(fifth and ninth columns of Table VII) because they cor- 
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respond to zone boundary modes in the (3 phase (see 
Table V) or because they are inactive T-point modes 
(see Tabic IV) agrees well with the results reported in 
Ref. [13]. 

The first anomaly is related to the experimentally ob- 
served IR mode with a frequency of 798 cm -1 in the a 
phase that does not disappear upon transition to the (3 
phase. Finnie et al. [44] explained this by suggesting that 
a two-phonon process in the (3 phase replaces the fun- 
damental mode of the a phase. Our calculations iden- 
tify two almost-degenerate IR modes that are close to 
this frequency, an A 2 mode at 751 cm -1 and an E mode 
752 cm -1 . In the a structure both of these modes are 
IR active, but the Born charge of the E mode is only 
0.02 while that of A 2 is 0.52, which means that the E 
mode in the a structure is almost invisible. In the f3 
structure, the A 2 mode disappears since it is no longer 
at r. On the other hand, the E mode remains at the 
r point and its Born charge is increased to 0.51. These 
results suggest a possible explanation for the "anomaly," 
namely that there are two IR modes in the a phase; one 
of them is much weaker than the other, but upon the 
transition to the (3 phase, the stronger one disappears by 
symmetry while the weaker one greatly increases its IR 
activity. The reason why the 752 cm -1 E phonon in the a 
structure acquires a larger Born charge upon converting 
to the (3 structure is that it gets some admixture of the 
440 cm -1 E mode, which has a much larger Born charge 
(0.74) than that of the 752cm -1 mode (0.02). 

The second anomaly is related to the 1076 cm -1 mode 
that is still present upon the transition to the (3 phase in 
the form of a fairly broad feature (see Fig. 1 in Ref. [43]), 
whereas it would be expected to vanish by symmetry. 
Swainson et al. [13] attributed this mode to a possible 
higher-order process. We think that it could also be re- 
lated to the fact that the 1050 cm -1 A 2 mode that is 
Raman inactive in the a-phase becomes Raman active in 
/3-phase. 

We also predict several phonon modes that are not seen 
experimentally, such as the Raman-active modes at 127, 
358, 378, 584, 752 and 1048 cm -1 . Since we have not 
computed Raman matrix elements, it is possible that the 
Raman intensities are small for these modes. We also 
find one weak IR-active mode at 752 cm -1 that is not 
seen in the experiments. 

A very low- frequency phonon at 50 cm -1 has been re- 
ported in a cristobalite. 42,45 ' 46 Wc believe this most likely 
corresponds to the Bi phonon that we have calculated to 
appear at 29 cm -1 , corresponding closely to the RUM 
mode that takes the ideal cubic cristobalite structure 
into the (3 structure. The same conclusion regarding 
the lowest-frequency Bi phonon was reached in Ref. [13]. 
The minimal-energy path between a and (3 phases that 
is related to this low-frequency phonon is discussed in 
Sec. IIIC2. 



C. SiC>2 cristobalite stability analysis 

As shown in Sec. IIIB, all calculated optical phonons 
in a and (3 cristobalites have lo 2 > 0, indicating that the 
relaxed structure is stable with respect to those modes. 
In view of the suggestion in Ref. [6] that the a phase 
might locally have D 2 rather than D4 point-group sym- 
metry, we checked carefully for instabilities leading from 
the a structure to P 2 structures, but found none. The 
possible subgroups of P4i2i2 (Df) having Z? 2 symme- 
try (without reduced translational symmetry) are C222i 
(Z? 2 ) and P2i2i2i (£>!), an d the phonon distortions lead- 
ing to these symmetry-lowered structures are the ones of 
P 2 and B\ symmetry respectively The lowest-frequency 
mode of P 2 symmetry is at 275 cm -1 , so there is certainly 
no sign of an instability there. On the other hand, the 
lowest-frequency B\ phonon is nearly soft at 29 cm -1 , 
suggesting that the a structure is nearly unstable to a 
spontaneous transformation into the P2i2i2i structure. 
To check this possibility more carefully, we started from 
the relaxed a structure and followed the distortion cor- 
responding to the 29 cm -1 B\ phonon, and confirmed 
that the energy increases monotonically (no double-well 
structure). Moreover, starting from one of these struc- 
tures having a small amount of the 29 cm -1 mode frozen 
in, a subsequent relaxation inside the resulting space 
group P2i2i2i lead to a recovery of the starting space 
group P4i2i2. We thus conclude, at least within our 
zero-temperature first-principles calculations, that the a 
structure is locally stable, i.e., does not spontaneously 
lower its symmetry to P> 2 . 

Nevertheless, the presence of several modes of quite 
low frequency in the a-cristobalite structure may be sug- 
gestive of low-energy pathways leading from the a to the 
(3 structure or between domains of the a structure. For 
example, we have shown above that the lowest-frequency 
29 cm -1 mode in the a-cristobalite structure corresponds 
to a phonon of the "ideal structure" that leads to the (3 
structure, and vice versa. This might suggest that there 
is a relatively low energy barrier in the configuration 
space that connects one structure to the other. Other 
phonons from the unstable triplet and doublet in the 
"ideal structure" have frequencies that arc higher, but 
still low enough to suggest that there is a low energy bar- 
rier for creation of the domains. In order to clarify these 
issues, we shall explore the energy landscape around the 
a and (3 structures in more detail in Sec. IIIC2. First, 
however, we begin with a general discussion of RUMs in 
the cristobalite phases in the next subsection. 

1. Rigid unit mode analysis 

Here we analyze the RUMs present in the high- 
symmetry cubic structure, 47 but constrained to maintain 
the periodicity of the four-formula-unit (-2=4) cell of the 
a-cristobalite structure (i.e., containing two unit cells of 
the /3-cristobalite structure). When these constraints are 



FIG. 3: (Color online) Five linearized RUMs in ideal cristobalite Z=4 structure. 



imposed, we find that there are five linearly independent 
RUM distortions as shown in Fig. 3. The first three of 
these distortions, which we denote as /3 2 , and (3% , con- 
sist of tetrahedral rotations of alternating signs about a 
single Cartesian axis, and carry the system into the I42d 
symmetry of the (3 phase. The last two, which we denote 
as «i and ol' 1i consist of a pattern of rotations around 
the [110] and [110] axes in the frame of Fig. 1, together 
with small translations of the tetrahedra needed to keep 
them connected at their apices, and carry the system into 
the P4i2i2 (or, for a[, into the enantiomorphic P4a2i2) 
symmetry of the a-cristobalite structure. Not shown in 
Fig. 1 are RUM rotations a.2 and a 2 associated with a 
second X point, and 013 and a' 3 associated with a third 
X point. 

Within the context of an ideal rigid-unit geometry (in 
which no additional relaxations are allowed), one can 
make the following mathematical analysis. The freezing 
in of the ideal (3\ RUM leads to a /3 structure oriented 
as in Sec. Ill A, whereas the freezing in of the ideal a\ 
RUM leads to an a structure as in that section. In the 
ideal $ structure, all five of the modes shown in Fig. 3 
remain as true RUMs - i.e., the tetrahedra can remain 
undistorted to first order in the mode amplitudes. Thus, 
all five modes arc expected to have low frequencies in a 
more realistic description. However, modes a.2, Q^j ®3i 
and a' 3 are no longer RUMs when a finite (3\ RUM is 
present. 

In the Si structure, only the &i, 5^, and f3\ distortions 
remain as true RUMs. However, the @2 and (3% modes at 
least share the same translational symmetry, and so may 
be expected to have somewhat low frequencies. The re- 
maining a2 , 0*2 > ^3 j an d modes are incompatible both 
in the RUM sense and in their translational periodicity. 48 

Not surprisingly, when we impose the translational pe- 
riodicity consistent with the five modes shown in Fig. 1, 
we confirm that these five distortions correspond quite 
closely to the five unstable phonon modes that we found 
in the ideal structure. The unstable T modes correspond 
to Pi , P2 and P3 , while the unstable X modes correspond 
to «! and a[. They also correspond closely to the low- 
frequency phonons in the a and f3 structures as discussed 
in Sec. IIIB3. 

Extending the mathematical analysis of the compati- 



bility of RUMs, it can be shown that there is an entire 
three-dimensional subspace of rigid- unit structures (i.e., 
with the tetrahedral rigidity condition satisfied exactly) 
in which finite rotations of type (ai,a4,/3i) are simul- 
taneously present, and having the space group P2i2i2i 
that is induced if any two of them are present. In a 
similar way, there are two additional three-dimensional 
(3D) manifolds (02, a'2^2) and (a3,a' 3 ,Ps) corresponding 
to different choices of the X point and thus having dif- 
ferent Z=A supercells. The three subspaces meet only 
at a single point (the cubic phase with all angles van- 
ishing), and RUMs selected from different 3D manifolds 
are always incompatible with each other in the sense that 
the perfect tetrahedral rigidity cannot be preserved when 
imposing both. This picture has important consequences 
for our understanding of the possible paths connecting 
domains of the a and (3 structures, as discussed below. 



2. Energy landscape inside 3D manifolds 

After we have explained the origin of the low-energy 
phonons in the a and (3 structures by relating them to 
RUM modes, we would now like to explore the energy 
landscape around these structures. To do so, we begin by 
finding a configuration space containing both structures. 
Since there is no group-subgroup relation between the a 
and (3 structures, we seek a maximal common subgroup 
of both structures. In the present case, this leads to the 
space group P2i2i2i (D^)- 

In the P2i2i2i configuration space, the a and (3 struc- 
tures represent two special points, and we know that the 
energy has local minima at these points because all com- 
puted phonon frequencies were found to be positive there. 
But then we also expect that there must be at least one 
saddle point connecting these points. To search for this 
saddle point, we started from the midpoint between the 
a and (3 structures in the 12-dimensional P2\l\l\ con- 
figuration space (described by nine internal coordinates 
and three cell parameters), and identified the unit vector 
e pointing between the two structures. We then applied 
a simple saddle-point search strategy in which compo- 
nent of the force vector parallel to e was reversed in sign 
before executing the steepest-descent update. This algo- 




FIG. 4: Energy in plane denned by the a structure (filled 
circle), $ structure (filled square), and saddle point (cross). 
Coordinates are chosen such that these structures occur at 
(0,0), (1,0), and (0.5,1), respectively. The energy difference 
separating contours is 3 meV per Si02 formula unit. 



rithm can be expected to succeed if the saddle point is 
not too far from the midpoint and if the principal axis 
of the negative Hessian eigenvalue at the saddle point 
is roughly parallel to e. In the present case, it led to 
a rapid convergence on the desired saddle point. Sur- 
prisingly, we find that the saddle point has a very low 
energy, only 5 meV per formula unit above that of the a 
structure, or 17 meV above that of the $ structure. 

The three points representing the a and (3 struc- 
tures and the saddle point determine a plane in the 12- 
dimensional configuration space. To confirm that the 
path running through the saddle point encounters only 
a single barrier, we have plotted the structural energy 
in this plane (without relaxation of other coordinates) in 
Fig. 4. We have somewhat arbitrarily carried out a lin- 
ear transformation on the coordinates in such a way that 
the a and (3 structures lie at (0, 0) and (1, 0) respectively, 
while the saddle point lies at (0.5, 1), in Fig. 4. The re- 
sults confirm the picture of a simple barrier of 5 meV 
encountered when going from the a to the f3 structure. 

Note that a transformation that would lead from the 
a to the f3 structure along a straight line in configura- 
tion space would have an enormously higher barrier of 
195 meV per formula unit. This is because the straight- 
line path is a poor approximation to a RUM. If instead we 
follow a curved minimum-energy path from a through the 
saddle to (3 and compute the relaxed Si-0 bond lengths 
and O-Si-0 bond angles along this path, we find that 
these remain almost constant. This strongly suggests 
that this minimum-energy path may be well approxi- 
mated by some RUM-like distortion. 

In the Sec. Ill CI, we pointed out that within the 
framework of ideal rigid-unit rotations, there is an en- 
tire three-dimensional subspace of structures for which 
the tetrahedral rigidity conditions are satisfied exactly, 
in which finite rotations of all three types arc present. 
We label an arbitrary configuration in this 3D manifold 
by (a\,a! 1 ,Px), where by convention the order of opcra- 
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FIG. 5: Energy as a function of rotation angles cj> ai and ^>p 1 , 
corresponding to rotations 5t\ and /3i shown in Fig. 3. The 
origin corresponds to "ideal" cristobalite. Filled squares at 
top and bottom denote /? minima, filled circles at left and 
right denote d\ minima, and crosses denote saddle points, as 
in Fig. 4. The energy difference separating contours is 3 meV 
per Si02 formula unit. 

tions is Qi followed by a! x and then /3i. The space group 
at a generic point in this 3D manifold is P2i2i2i, the 
same one we have just been discussing. It thus seems 
likely that the minimum-energy path in Fig. 4 may cor- 
respond approximately to a path from the point (ai,0,0) 
to the point (0,0, /3i) and lying, at least approximately, 
in the two-dimensional (2D) subspace (ai,0,,3i). 

To test this conjecture, we first created an ideal rigid- 
unit structure for each pair of angles (ai,/3i) on a two- 
dimensional mesh. We then used our first-principles cal- 
culations to relax each structure subject to the constraint 
that these two angle variables should not change. Tech- 
nically, we did this by carrying out the minimization of 
the energy in the ten-dimensional subspace orthogonal 
to the two-dimensional surface for each starting point 
(ai,/3i). We typically found that these relaxations were 
small, confirming the approximate validity of the RUM 
picture. 

The energy surface determined in this way is plotted 
as a function of rotation angles a± and /3i in Fig. 5. 
The minima corresponding to the a structure are im- 
mediately visible near the left and right sides of the 
figure, while those corresponding to the slightly lower- 
energy (3 structure appear near the top and bottom. The 
minimum-energy path appears to be roughly circular on 
this plot, and four equivalent saddle points are appar- 
ent at ai ~ ±19° and j3\ ~ ±13°. These saddle points 
are equivalent to the one identified in Fig. 4, with a bar- 
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rier height of 5 meV per formula unit relative to the a 
structure. We thus confirm the presence of a very low- 
energy barrier between these structures, and identify it 
as approximating a certain path in the space of rigid- 
unit rotations. A video animation showing the evolution 
of the structure along this path is provided in the sup- 
plementary material. 49 

It is important to note that, according to the simplified 
model of Eq. (1) of Ref. [7], the energy would remain ex- 
actly zero on the entire («i, surface of Fig. 5 since the 
ideal rigid-unit structures satisfy the rigidity conditions 
analytically. The RUM framework envisages extensions 
to make the model more realistic; one way to do this is by 
adding an energy term that depends on the relative tilts 
of neighboring tetrahedra. 50 We tried this by introduc- 
ing a simple double-well potential model that penalizes 
the departure of the Si-O-Si bond angles from a preferred 
bending angle. In this model the change of total energy 
per formula unit is 



AE = 




where the sum runs over all N Si-O-Si bond angles 4>i 
in the unit cell. We found that we could obtain an opti- 
mal fit 51 to the results of our first-principles calculations 
using parameters Eq = 83 meV per formula unit and 
4>q = 145°. The energy landscape of the fitted model 
looks very similar to the results plotted in Fig. 5. In par- 
ticular, the overall circular aspects of the energy land- 
scape and minimal-energy path in Fig. 5 are reproduced. 
However, the model unfortunately assigns identical ener- 
gies to the <5i and (3\ structures, and moreover predicts a 
path connecting them on which the energy remains com- 
pletely flat. This happens because, for any given pair of 
angles on or near this path, one can find a small 

cti such that the rigid-unit structure (ai,a[,f3i) has all its 
Si-O-Si bond angles exactly equal to 4>q. Therefore our 
simplified model of Eq. (7), or any other model that de- 
pends solely on the Si-O-Si angles, predicts a zero-barrier 
path between ot\ and j3\ structures. This behavior is rem- 
iniscent of an early model of Nieuwenkamp 52 for (3 cristo- 
balitc, in which the Si-O-Si bond was assumed to rotate 
freely on an annulus lying in the plane that is equidistant 
between Si atoms. 

To further test the model of Eq. (7) , we performed first- 
principles calculations on a mesh of (unrelaxed) struc- 
tures in the 2D space (ai,a' 1 ,/?i=0). The resulting energy 
landscape is shown in Fig. 6. The barrier between a\ and 
a[ structures is now about 35 meV, substantially higher 
than for the path connecting a and (3 structures. When 
we use the same fitting parameters obtained earlier, we 
again get very good overall agreement; the energy land- 
scape obtained from our model has the same diamond- 
like appearance as in Fig. 6, and saddle points appear in 
very similar locations. Moreover, the barrier of 41meV 
predicted by the model is in quite good agreement with 
the first-principles value of 35 meV. However, in this case 
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FIG. 6: Energy as a function of rotation angles <j> ai and <j> a i , 
corresponding to rotations <Si and shown in Fig. 3. The 
origin corresponds to "ideal" cristobalite. The four minima 
(filled circles) correspond to the a[ structure (top and bot- 
tom) and to the &i structure (left and right). The energy 
difference separating contours is 3 meV per Si02 formula unit. 

the picture presented by Fig. 6 is somewhat misleading, 
because it turns out that the entire minimum-energy path 
lying in the /?i=0 plane is unstable, and falls to lower 
energy as (3\ is turned on. Thus, the apparent saddle 
points in Fig. 6 arc actually stationary points with two 
negative eigenvalues in the 3D (ai,a[,[3i) space. Within 
the model of Eq. (7), in fact, the lowest-energy path con- 
necting the a.\ and a' x structures is actually completely 
flat, being composed of a segment connecting a\ to (3\ 
and then another connecting (3\ to a^. This observation 
agrees with our first-principles calculations, since if we 
start from the purported saddle-point configuration and 
do a structural relaxation subject to the constraint that 
ai = a' l7 the structure is found to converge to the /?i 
structure as expected. 



3. Cell volume at minima and saddle point 

Because we have found the unit-cell parameters to be 
very sensitive to details of the calculation, we increased 
the energy cutoff from 22 Ha to 30 Ha in order to obtain 
an accurate description of the volume changes along the 
minimum-energy path. We obtain a volume per formula 
unit of 45.7 and 46.7 A 3 for the a and /3 local minima re- 
spectively, so that the volume is about 2.2% larger for the 
latter. This is in qualitative agreement with experiments, 
which show that the (3 structure is about 5% larger 53 , and 
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implies that applied pressure would tend to favor the a 
phase and raise the a-to-/3 transition temperature. At 
the saddle point, we find that the volume per formula 
unit is 46.8 A 3 , which is just slightly larger than for ei- 
ther of the parent-phase structures. This finding may be 
of interest for future studies of the pressure-dependence 
of the phase-transition mechanism. 



4- Domain walls 

The barriers discussed in Sec. Ill C 2 refer to transfor- 
mation pathways in which the crystal remains periodic 
and transforms homogeneously, and the energy barriers 
are given per unit cell. It would also be of interest to con- 
sider the energies of domain walls between various a and 
(3 structures. This is beyond the scope of the present in- 
vestigation, but the results for homogeneous transforma- 
tions may give some hints as to what could be expected. 
For example, we speculate that domain walls connecting 
a and (3 structures belonging to the same 3D rigid-unit 
manifold will probably have a rather low energy per unit 
area, while those connecting structures belonging to dif- 
ferent 3D manifolds would be expected to have much 
higher energies. 



IV. DISCUSSION 

In this section we give a brief overview of several pre- 
viously proposed models of a and (3 cristobalite phases, 
and discuss how the results of our calculations relate to 
those models. 

The RUM model of Ref. [7] describes the (3 phase as an 
average cubic structure that has strong dynamical fluctu- 
ations occurring simultaneously into RUMs in all allowed 
regions of the Brillouin zone. A simplified version of this 
picture would be one in which the tetrahedra are assumed 
to be completely free to pivot around their apices, as in 
Eq. (1) of Ref. [7]. In general, the simultaneous excita- 
tion of more than one RUM will have an associated en- 
ergy cost because the tetrahedra typically cannot remain 
perfectly rigid while undergoing both kinds of distortion 
simultaneously. However, as an exception, we have iden- 
tified 3D rigid-unit manifolds within which the geomet- 
rical constraints can simultaneously be satisfied. Within 
the model of Eq. (1) of Ref. [7], or the split-atom model 
of Ref. [54] , the energy landscape within this special 3D 
manifold would be completely flat, and one would expect 
that freezing in of one RUM of type fii, a[ or (3\ would 
have no consequence on the energy profile of one of these 
other RUM distortions. 

However, once one goes beyond the simplest versions of 
the model and includes terms that depend on the Si-O-Si 
bond angles at the apices, our calculations indicate that 
the RUM distortions of type a.\, a[ and f3\ become cou- 
pled and have a rich energy landscape. As a step in this 
direction, the more sophisticated split-atom model hav- 



ing an additional energy term depending on Si-O-Si bond 
angles 50 should provide an improved description. How- 
ever, we note that even this model, or any model based 
solely on Si-O-Si bond angles, still has a nonphysical be- 
havior in that it would necessarily predict zero-energy 
barriers between the a,\ and (3\ structures, as discussed 
at the end of Sec. IIIC2. Nevertheless, we believe that 
the split-atom and similar models can provide important 
complementary information to ours, since they are not 
restricted to periodic supercell structures as ours are. 

Among the models of cristobalite phase transitions 
is also the model of Hatch and Ghose. 55 They argue 
that the (3 phase is dynamically and spatially fluctuat- 
ing between the 12 different possible a domains having 
P4i2i2 space-group symmetry. The counting arises be- 
cause there are three different X points; each exhibits a 
doublet of degenerate modes leading to enantiomorphic 
a\ and a[ structures (see Fig. 3); and the tetrahedra 
can rotate by ±4>. The model is based on symmetry 
arguments and assumes that all of the barriers separat- 
ing these 12 a structures are small. However, our work 
suggests that the barriers separating different types of 
a domains have very different barriers. Furthermore, 
their model does not take into account the fact that the 
(3 structure is easily accessible with a very low barrier, 
suggesting that fluctuations into the (3 structure may be 
more important than some of the other a structures. 

Finally, O'Keeffe and Hyde 4 do discuss a path connect- 
ing a and (3 structures, but it is of a different type than 
those discussed above since it connects a and (3 struc- 
tures belonging to different 3D rigid-unit manifolds. In 
our notation, their path would connect a± or a[ to fii 
or /?3, etc. Such a path would involve the simultaneous 
application of RUM rotations that are incompatible with 
each other, and as such would be expected to have a high 
energy barrier. 

To clarify our view of the cristobalite phase transitions, 
we start by emphasizing once again the existence of three 
distinct 3D rigid-unit manifolds, as described above at 
the end of Sec. Ill CI. To review, one of these is de- 
scribed by rotation angles (ai,a[,f3i) giving rise to struc- 
tures of space group P2i2i2i whose translational period- 
icity is that corresponding to the X point (27r/a c )(001) 
or cquivalently (27r/a c )(110). This manifold contains 
the ai, a[, and (3\ structures, and their partners with 
reversed sense of rotation, as shown schematically in 
Fig. 7. The second and third 3D subspaces are de- 
scribed by rotations {a.2,a.' 2 ,f32) and ((23,0:3, /3s), with pe- 
riodicities set by X points (27r/a c )(010) = (27r/a c )(101) 
and (27r/a c )(100) = (27r/a c )(011), respectively We have 
found that these three subspaces are essentially incom- 
patible, in the sense that it is not possible to combine 
rotations taken from any two of them into a combination 
that preserves the rigid-unit constraints. This occurs in 
part because these three 3D rigid-unit subspaces have in- 
compatible translational symmetries, but also because of 
incompatibilities in the patterns of rotations. 

The structure of the space sketched in Fig. 7 is intended 
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FIG. 7: Sketch of important states in one of the three- 
dimensional rigid-unit subspaces discussed in the text. Local 
energy minima associated with a (D4) and (3 (D2d) struc- 
tures are indicated by filled circles and squares respectively. 
Remainder of space, including saddle points (crosses), has D2 
symmetry. 

to reflect a three-level hierarchy of energies and energy 
barriers as suggested by our analysis. In the model of 
Eq. (7), the energy is degenerate for all six of the struc- 
tures shown in Fig. 7, as well as on the solid curves con- 
necting them. According to our first-principles results, 
this picture is modified so that the a and (3 structures are 
local minima, with low-energy saddle points (~5 meV) 
between them (see Fig. 5). The low curvature of the en- 
ergy surface along these curves is reflected in the presence 
in Table VI of a very soft 29 cm -1 Bi mode starting from 
the a structure, and a 35 cm^ 1 M3M4 doublet starting 
from the structure. 56 While our calculations have the 
a structures at a slightly higher energy than the (3 ones, 
this is presumably reversed in the true physical system. 

The next energy scale in the hierarchy is that asso- 
ciated with the direct paths between a structures in 
the same 3D manifold, indicated by the dashed lines in 
Fig. 7. As shown in Fig. 6, this energy is on the order 
of ^35 meV, so that the true minimum-energy path be- 
tween neighboring a structures goes instead through (or 
perhaps nearly through) the (3 structures. 

Finally, the highest energies are associated with the 
barriers separating any of the structures in Fig. 7 from 
any of the structures in the other two 3D subspaces. 
These barriers are on the order of 80meV, the energy 
needed to pass through the undistorted cubic phase. 
While not enormously larger than the 35 meV mentioned 
above, this is high enough that we do not expect these 
barriers to be especially relevant for the phase transitions 
in this system. 

We can now speculate on the nature of the phase tran- 
sition between a and j3 cristobalites. We propose that 
in the lower-temperature a phase, the system is locally 



frozen onto one of the minima of type a in one of the 
3D manifolds, but with substantial fluctuations along 
the low-energy paths leading to the two neighboring (3 
structures in the same manifold. Then, in the higher- 
temperature (3 phase, we speculate that the system in- 
stead shifts over and condenses locally onto one of these 
(3 structures, but with substantial fluctuations along the 
low-energy paths leading to the four neighboring a struc- 
tures, all in the same 3D manifold. The fact that there 
are four low-energy paths to fluctuate along, instead of 
two, is consistent with the fact that the (3 phase (be- 
ing the higher-temperature phase) has higher entropy. If 
the system were truly to freeze onto a single (3 struc- 
ture, it would be globally tetragonal, with space group 
IA2d. However, it is also possible that the system forms 
on some larger scale into spatiotemporal domains com- 
posed of (3 structures from all three of the 3D manifolds, 
giving an overall average Fd3m structure in accord with 
the picture espoused in Refs. [7,8]. 

Let us return for a moment to the recent work of Zhang 
and Scott, 6 who argued that their Raman studies of sin- 
gle crystals of /3-cristobalite were inconsistent with Oh 
symmetry. Assuming Did symmetry instead for the f3 
phase, these authors then noted that D4 is not a sub- 
group of Did-, and thus that the existence of a group- 
subgroup relation for the phase transition would rule out 
the assignment of the a phase to the D4 a structure. On 
this basis, they suggested that a lower symmetry, such 
as D 2 , should be considered for a-cristobalite. Our view, 
instead, is that a group-subgroup relation does not have 
to hold for the transition, since the transition is known 
to be of first order, and thus assignments of Did and 
Z?4 for the a and (3 phases respectively are not inconsis- 
tent. As pointed out in the introduction, while certain 
spectroscopic signatures of the transition are indicative 
of a weakly first-order transition, the volume change and 
latent heat at the transition are substantial. The transi- 
tion may perhaps be described as a reconstructive transi- 
tion in the sense of Tolcdano and Dmitriev, 11 although in 
the present case the rearrangements of atoms can occur 
very gently, because of the existence of very low-barrier 
paths of Z?2 symmetry connecting the D4 (a) and Did (/?) 
structures. The situation may be somewhat analogous 
to the tetragonal-to-orthorhombic and orthorhombic-to- 
rhombohcdral transitions in ferroelectric perovskites such 
as BaTiOs and KNb03, where the presence of low-barrier 
paths of monoclinic symmetry is associated with the 
weakly first-order nature of the transitions. 57 

Unfortunately our calculations are carried out at K 
with crystal periodicity imposed. It is therefore diffi- 
cult to draw any firm conclusions about the nature of 
the phase transitions between cristobalite phases, espe- 
cially if fluctuations are as important as we think they 
are, and much of what we have said above must remain 
speculative. Nevertheless we hope that the results of our 
calculations will be of use in developing improved models 
that may allow for realistic unite-temperature modeling 
of the phase transitions in this system, ultimately leading 
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to a resolution of the controversies that have surrounded 
this system over the years. 



V. SUMMARY 

Based on first-principles calculations, we have per- 
formed a detailed analysis of the a (P4i2i2) and (3 (I42d) 
structures of cristobalitc SiC>2- In particular, we have 
confirmed that both structures are locally stable against 
all possible distortions associated with r-point modes of 
the four-formula-unit conventional cell. We have calcu- 
lated phonon frequencies for the a and (3 structures, com- 
pared these to the experimental values, and discussed 
how the phonons in these two structures are related to 
each other. We have also tried to resolve some experi- 
mental anomalies that were found in spectroscopic stud- 
ies of the cristobalite phases. Finally, we have explored 
the energy landscape connecting the a and (3 structures. 
We have emphasized the existence of three distinct 3D 
manifolds of structures, each of which contains both a 
and (3 structures that can be connected to each other 
within the manifold by paths with a surprisingly small 
barrier of 5 meV per formula unit, while paths connecting 
different manifolds have a much higher barrier. While our 
calculations do not properly treat fluctuations, we nev- 
ertheless have speculated on the possible consequences 
of our findings for the understanding of the a-(3 phase 
transition in cristobalite SiC>2. 
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APPENDIX A: COMPARISON WITH 
RIGID-UNIT GEOMETRY 

In a picture in which the rigid-unit constraints are per- 
fectly enforced, it turns out that the structures of a and 
(3 symmetries are completely determined by a single pa- 
rameter, which can be taken to be the volume V per for- 
mula unit relative to the corresponding value Vq in the 
ideal cubic structure. (That is, Vq is the volume below 
which rigid distortions start to appear, as explained in 
Sec. Ill A.) In this Appendix, we check to see how closely 
our structures, as optimized from the first-principles cal- 
culations, match with this elementary model. 

The solid curves in Figs. 8 and 9 show the mathe- 
matical predictions of this elementary model, obtained 
by applying rotations of types ce\ and (3\ (see Fig. 3) 
in such a way as to keep the tetrahedra perfectly rigid. 
(For V > Vq, the elementary model cannot be satisfied, 
and the ideal cubic parameters arc plotted instead.) The 
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FIG. 8: Structural parameters of the a structure vs. vol- 
ume per formula unit. Top panel: c/a ratio. Bottom panel: 
absolute values of deviations of internal parameters u(Si) 
(squares), x(0) (triangles), y(0) (crosses), and z(0) (circles) 
from ideal-cubic values. Symbols represent first-principles cal- 
culations; lines are fits to an ideal rigid-unit geometry. 



symbols shown in Figs. 8 and 9 denote the results of 
our first-principles calculations where, for each specified 
value of V, the volume was treated as a constraint while 
all other structural parameters were relaxed. The fit was 
optimized by choosing a common V = 55.1 A 3 for both 
a and (3 structures. For reference, the first-principles 
equilibrium volumes are 45.7 and 46.7 A 3 for the a and 
(3 structures, respectively. 

We find that the agreement is extraordinarily good for 
all of the internal parameters, but that there are some sig- 
nificant discrepancies in the c/a ratios. At first sight this 
may seem contradictory: why are the c/a ratios off by 
many percent, while the Si-0 bond lengths agree within 
^0.05%? The answer is connected with the presence of 
volume-preserving tetragonal distortions of low energy 
cost. In such a distortion, each tetrahedron is stretched 
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FIG. 9: Structural parameters of the structure vs. volume 
per formula unit. Top panel: c/a ratio. Bottom panel: x(0). 
Symbols represent first-principles calculations; lines are fits 
to an ideal rigid-unit geometry. 
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slightly along c and compressed in a (or vice versa), and 
it happens that the tetrahedral angle of arccos(l/v / 3) is 
precisely the one at which Si-0 bond lengths are pre- 
served to hrst order in the distortion amplitude. While 
the O-Si-0 bond angles do change at first order, this 
may entail a smaller energy cost than for bond-length 
changes. As expected from this analysis, we find that our 



first-principles O-Si-0 bond angles differ from the ideal 
ones by ~4%. In short, it appears that it is energetically 
more important to preserve bond lengths than bond an- 
gles, and that for geometrical reasons this translates into 
an enhanced freedom for the c/a ratio to deviate from 
the ideal rigid-unit geometry. 
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